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Abstract 

The causal interpretation of quantum mechanics is applied to a homogeneous and isotropic 
quantum universe, whose matter content is composed by non interacting dust and radiation. For 
wave functions which are eigenstates of the total dust mass operator, we find some bouncing 
quantum universes which reachs the classical limit for scale factors much larger than its minimum 
size. However these wave functions do not have unitary evolution. For wave functions which are 
not eigenstates of the dust total mass operator but do have unitary evolution, we show that, for 
flat spatial sections, matter can be created as a quantum effect in such a way that the universe can 
undergo a transition from an exotic matter dominated era to a matter dominated one. 
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I. INTRODUCTION 



The Bohm-de Broglie (BdB) interpretation [l| has been sucessfully applied toquan- 
tum minisuperspace models llo| . and to full superspace |rj] |l3j]. In 



the first case, it was discussed the classical limit, the singularity problem, the cosmo logical 

constant problem, and the time issue. It was shown in scalar field and radiation models 

for the matter content of the early universe that quantum effects driven by the quantum 

potential can avoid the formation of a singularity through a repulsive quantum force that 

counteract the gravitational attraction. The quantum universe usually reach the classical 

limit for large scale factors. However, it is possible to have small classical universes and 

n 

large quantum ones: it depends on the state vector and on initial conditions |9|. It was 
also shown that the quantum evolution of homogeneous hypersurfaces form the same four- 
geometry independently on the choice of the lapse function 

In the present work we study the minisuperspace model given by a quantum Friedmann- 
Lemaitre-Robertson- Walker (FLRW) universe filled with dust and radiation decoupled from 
each other. We write down the hamiltonian that comes from the velocity potential Schutz 
formalism Q|. After implementing a canonical transformation, the momentum associated 
to the radiation fluid pr and to the dust fluid p v appear linearly in the superhamiltonian 
constraint. Both can be associated to time parameters, but physical reasons and mathemat- 
ical simplicity led us to choose the coordinate T associated with p? as the time parameter. 
This is equivalent to choose the (reversed) conformal time. We quantize this system obtain- 
ing a Schrodinger-like equation. We analyze its time dependent solutions applying the BdB 
interpretation in order to study the scale factor quantum dynamics. 

We first consider an initial quantum state given by a gaussian superposition of the scale 
factor which is an eigenstate of the total dust mass operator (matter is not being created 
nor destroyed in such states), and we compute the solution at a general subsequent time 
by means of the propagator approach. We calculate the bohmian trajectories for the scale 
factor. For flat and negative curvature spatial sections, we find that the quantum solutions 
for the scale factor reach the classical behaviour for long times, but do not present any initial 
singularities due to quantum effects. In the same way, in the case of positive curvature spatial 
sections, the classical initial and final singularities are removed due to quantum effects, and 
the scale factor oscillates between a minimum and a maximum size. For large scale factor, 
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the classical behaviour is recovered. However, such eigenf unctions of the total dust mass 
operator do not have unitary evolution. This led us to consider an initial state given by 
gaussian superpositions of the total dust matter content. In this situation, dust and radiation 
can be created and destroyed. We calculate general solutions for flat, negative and positive 
curvature spatial sections. In particular, for flat spatial sections, we construct a wave packet 
whose quantum trajectories represent universes which begin classically in an epoch where 
the dust matter has negative energy density (exotic dust matter), evolving unitarily to a 
configuration where quantum effects avoid the subsequent classical big crunch singularity, 
performing a graceful exit to an expanding classical model filled with conventional matter 
and radiation. There is thus a transition from an exotic matter era to a conventional matter 
one due to quantum effects. 

This paper is organized as follows. In section |H] we synthesize the basic features of the 
Bohm-de Broglie interpretation of quantum mechanics, which will be necessary to interpret 
our quantum model studied in other sections. In section IIH[ we briefly summarize the 
velocity potential Schutz formalism, and we apply it to construct the hamiltonian of the 
FLRW universe filled with two perfect fluids, which are dust and radiation. We then review 
and analyze the classical features of the two perfect fluids FLRW model in order to have 
the results to be contrasted with the quantum models of the following sections. In section 
IIV1 we present some new results concerning the existence of singularities in the quantization 
of the one fluid case. We show that, when the fluid is radiation, all quantum solutions 
do not present singularities. In section El we quantize the model with two fluids, and we 
compute the solutions of the Schrodinger like equation for two different initial conditions: 
the first being an eigenstate of the total dust matter operator, and the second a gaussian 
superposition of total dust matter eigenstates. We interpret the solutions according to the 
BdB view and we develop the main results of the paper. Section IVII is for discussion and 
conclusions. 



II. THE BOHM-DE BROGLIE INTERPRETATION OF QUANTUM MECHAN- 
ICS 

In this section, we briefly review the basic principles of the Bohm-de Broglie (BdB) 
interpretation of quantum mechanics. According to this causal interpretation, an individual 
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physical system comprises a wave ^>(x,t), which is a solution of the Schrodinger equation, 
and a point particle that follows a trajectory x(t), independent of observations, which is 
solution of the Bohm guidance equation 



p = mx = VS(x,t)\ x=x ( t) , (1) 

where S(x,t) is the phase of In order to solve Eq.(U), we have to specify the initial 
condition x(0_) = x . The uncertainty in the initial conditions define an ensemble of possible 
motions, 

It is sufficient for our purposes to analyze the Schrodinger equation for a non relativistic 
particle in a potential V(x), which, in coordinate representation, is 



dV(x,t) 
i- 



i v2 + ^» 



Substituting in (J2J) the wave function in polar form, \I/ = Aexp(iS), and separating into real 
and imaginary parts, we obtain the following two equations for the fields A and S 

dt 2m 2m A w 

Equation © can be interpreted as a Hamilton- Jacobi type equation for a particle submitted 
to a potential, which is given by the classical potential V(x) plus a quantum potential defined 
as 

It is possible to verify that the particle trajectory x(t) satisfies the equation of motion 

m $ = - W - VQ. (6) 

The classical limit is obtained when Q — 0. The BdB interpretation does not need a 
classical domain outside the quantized system to generate the physical facts out of poten- 
tialities. In a real measurement, we do not see superpositions of the pointer apparatus 
because the measurement interaction causes the wave function to split up into a set of non 
overlapping packets. The particle will enter in one of them, the rest being empty, and it will 
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be influenced by the unique quantum potential coming from the sole non zero wave function 
of this region. The particle cannot pass to another branch of the superposition because they 
are separated by regions where ^ = 0, nodal regions. 

In section |3 the FLRW minisuperspace model containing dust and radiation as two 
perfect decoupled fluids will be quantized. A preferred time variable can be chosen as one 
of the velocity potencials associated to the fluids (radiation), yielding a Schrodinger like 
equation. Then the BdB interpretation of our quantum model runs like in Ref . , in close 
analogy to the non relativistic particle model described above. In the present case, however, 
the scale factor of the universe will not be the only degree of freedom: the velocity potential 
associated with the dust field and its canonical momentum, interpreted as the dust total 
mass, are also present. They satisfy a Hamilton- Jacobi equation modified by an extra term, 
the quantum potential, so that their time evolution will be different from the classical one. 
The main features of this classical model we describe in the following section. 



III. CLASSICAL DUST PLUS RADIATION MODEL IN THE VELOCITY PO- 
TENTIAL SCHUTZ FORMALISM 

We start by considering a perfect fluid in a FLRW universe model. The line element is 
given by 

ds 2 = -N 2 dt 2 + a 2 (t) lij dx i dx j (7) 

where N is the lapse function, a is the scale factor, and 7^- is the metric of the three- 
dimensional homogeneous isotropic static spatial section of constant curvature k = 1,0, or 
-1. 

Following the Schutz's canonical formalism to describe the relativistic dynamics of a 
perfect fluid in interaction with the gravitational field Q|, we introduce the five velocity 
potentials, a, (3, 8, (p and s. The potentials a and (3, which describe vortex motion, vanish 
in the FLRW model because of its symmetry. The potential s is the specific entropy and 9 
can be related with the temperature of the fluid. By now (p works only as a mathematical 
tool. 

The four-velocity of the fluid is obtained from the velocity potentials as 
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U v = -{(p,„+6s, v ), (8) 
where [i stands for the specific enthalpy. The four velocity is normalized as 

g a pU a U p = -l. (9) 

Using this equation, it is possible to write the specific enthalpy ji as a function of the velocity 
potentials. 

The action for a relativistic perfect fluid and the gravitational field in the natural units 
c = h — 1 is given by 

I = -^nf d A x^ A K+ f d 4 x^g-p+^-[ d 3 xVhh i3 K*\ (10) 
ot* Jm Jm 3l£ JdM 

where l p = (SttG/S)^ 1 ^ 2 , G is the Newton's constant (hence l p is the Planck length in the 

natural units), A 1Z is the scalar curvature of the spacetime, p is the pressure of the fluid, 

hij is the three metric on the boundary dM of the 4-dimensional manifold M, and 

its extrinsic curvature. The velocity potentials are supposed to be functions of t only, in 

accordance with the homogeneity of spacetime. The perfect fluid follows the equation of 

state p = \p. 

Substituting the metric (JJJ) into the action (jlOj). using the formalism of Schutz ^| to 
write the pressure of the fluid as 



V = POr 



(p + 6s 



A+l 

A 



ex P(-7i:)> ( n ) 

OQ T A, 



_N{\ + 1) 

with p 0r and s 0r constants, computing the canonical momenta p v ,p s ,Pe for the fluid and 
p a for the gravitational field, using the two constraints equations pe = 0, 9p v = p s , and 
performing the canonical transformation 

r=-^e*p(--W A+1) PUn (12) 

and 



= <f + (A + l)s 0r — , (13) 

Pip 



leading to the momenta 



7 (A+1) 
jP 

POr ^ S 0r 



p^e^^V-expf^-), (14) 



and 

P VN =Pv, (15) 
we obtain for the final Hamiltonian (see Ref. [16] for details), 

^ m = Hi- fa+ S' (16) 

where N plays the role of a Lagrange multiplier whose variation yields the constraint equa- 
tion 

(17) 

where ~ means 'weakly zero' (this phase space function is constrained to be zero, but its 



Poisson bracket to other quantities is not). We have redefined a = yV/(lQirlp) a in order 
for a be dimensioless, and N = where V is the total comoving volume of the spatial 

sections. The tilda have been omitted. Considering now two decoupled fluids, one being 
radiation (A r = 1/3), and the other dust matter (A^ = 0), the Hamiltonian reads: 

H = NH = N^^-Qna + ^+p^j (18) 

The classical Hamilton equations are: 

r TT^ N Ylaa 
a = {a,H\ = -— p a ^p a = — , (19) 

P* = {P*,H} = n(-^ + 6k + %), (20) 
N 

T=-, (21) 
a 



V = &,H} = N, (22) 

Pt = Pip = =r- Pt, p v are constants. (23) 

The superhamiltonian is constrained to vanish due to variation of the Hamiltonian with 
respect to the lapse function N, Ji ~ 0, 

2 

''<• QKa + ^+ P(p = 0. (24) 



24a a 
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The constraint (ffifj) combined with Eqs. (|19p and (}2~3*)) yield the Friedmann's equation 



N 2 



a 2 6 \a 4 a 3 



(25) 



Note that the conjugate momenta pt and p^, classical constants of motion, can be identified 
to the total content of dust and radiation in the universe: 



p v = WnGa^pm, (26) 

p T = 16nGa 4 p r . (27) 

Note also that Eq. (j22j) implies that dtp = Ndt, hence p is cosmic time, while Eq. (j21j) 
yields adT = Ndt so T is conformal time. Consequently, choosing N = 1 means taking 
coordinate time t as cosmic time (p, while choosing N = a imposes coordinate time to be 
conformal time T. Explicit analytic solutions of Eqs. 1)1912012312 5ft can be obtained only in 
the gauge N = a. In this gauge, besides the constraint with N = a, we obtain the 
simple second order equation, 

(28) 



a" + KCL 



where a prime means differentiation with respect to conformal time, which we denote r] from 
now on. The solutions read: 



[1 - cosfa) +r] eq sm(ri)} 



cy 



2o < 



2^L + 


fjL\ 2 


Veq 





1, 



; k = 0, 



(29) 



-1. 



[cosh(?7) + r] eq smh(r]) — 1] ; k 

The quantity a eq is defined to be the value of the scale factor at the equilibrium time where 
p m = p r , and rf eq = 3/ (2tc G p r a A ) = 24a eq / \ p v |. 

As we will see in section (jVj), the presence of quantum effects can create exotic dust 
matter content. Hence, for comparison, we analyze a classical universe filled with exotic 
dust, which means p m < 0, i.e p v < 0. For simplicity, let us focus on the flat spatial case. 

In the presence of exotic dust, the behaviour of the scale factor is drastically different. 
From the Friedmann Eq.(J2SJ), since p^ < 0, the radiation density must always be equal or 
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greater then the dust density, otherwise the Friedmann equation 

2 



1 

6 I a 2 



has no solution. For small values of the scale factor, the radiation term dominates. As the 
scale factor grows, the exotic dust term begins to be comparable to the radiation term up 
to the critical point where both are equal and a = 0. ^From this point, the scale factor 
decreases until the universe recollapses. Note that Eq.(j2Hj) for k = and p v < implies 
that a" < at all times. Hence, contrary to the normal dust matter case where after the 
big bang the universe expands forever [see Eq. (|29j) for k — 0], in the exotic case the universe 
recollapses in a big crunch. The qualitative evolution of the scale factor is plotted in figure 
1: The deceleration parameter in conformal time is given by 




FIG. 1: Scale factor evolution of a universe filled with exotic dust (p v < 0). 



a"a 



+ 1. 



(30) 



It diverges when the scale factor reaches its maximum value (a' = and a" < 0). 



IV. FLRW QUANTUM MODEL WITH RADIATION 



In this section, we present a general result concerning the presence of singularities in the 
quantization of a FLRW model with radiation. The hamiltonian constraint in this case is 



and rj is conformal time, as discussed above. 

Using the Dirac quantization procedure, the hamiltonian constraint phase space function 
7i becomes an operator which must annihilate the quantum wave function: 7^\& = 0. One 
then obtains in natural units the Wheeler-De Witt equation for the minisuperspace FLRW 
metric with radiation: 

(a ' v)= {~Ta& + 6 ™ 2 ) * (a ' v) ■ (32) 

Note that a particular factor ordering has been chosen and, because p v appears linearly in 
Eq.(|31j). rj — > —77 is chosen to be the time label in which the wave function evolves (the sign 
rever sing was done in order to express this quantum equation in a familiar Schroedinger 
form |lfij|). 

The scale factor is defined only in the half line [0, oo), which means that the superhamil- 
tonian (}3"T]) is not in general hermitian. Hence, if one requires unitary evolution, the Hilbert 
space is restricted to functions in L 2 (0, oo) satisfying the condition 

— (0,^=0^(0,7?), (33) 

where a is a real parameter j^ . 

We will now show that condition ()33j) . together with the assumption that \]/(a, 77) is 
analytic in rj at a = 0, implies that general quantum solutions of Eq. (j32J) . when interpreted 
using the BdB interpretation, yield quantum cosmological models without any singularity. 

We can rearrange Eq.(j22J) in order to isolate the second spatial derivative: 



q2 

r * (a, 77) =24 



d 

-i— # (a, rj) + 6 k a 2 ^ (a, rj) 
or] 



(34) 



da 2 

Using the BdB interpretation, the scale factor equation of motion is given by the gradient 
of the phase S (a, rj) of the wave function 

a =-S a (a,r,) = -- — = / (a, rj) , 

where the index a means derivative with respect to a. Taking the boundary condition (133(1 
at a = 0, the velocity function / (0, 77) vanishes. Hence, if there is a time 770 where a(r] Q ) = 0, 
then a'(r]o) — 0. For a" one has: 

a » = ?L a > + ?L = ?L f + ?L. 

da dr] da drj 
10 



This is also zero at a = unless df /da diverges there. However, 



df 
da 



Van®* 



24 



+ 



(^^:) 2 - ( v 



24 (^^*) 2 



is obviously finite if condition (jHHj) and analyticity of ^ in rj is satisfied at a = 0. The case 
when (\I/\I/*) 2 = does not need to be analyzed because bohmian trajectories cannot pass 
through nodal regions of the wave function. 

The same reasoning can be used for all higher derivatives d n a/dr] n at a = to show that 
they are all zero: one just have to use equation ()34j) to substitute d 2 ^ /da 2 for d^//dr] and 
condition (|3*3|) to substitute d 2 ^ /dadrj for ad^//dr] at a = whenever they appear, and 
then use analyticity of \I/ in rj at a = 0. 

With these results, if there is a time t)q where a(r] Q ) = 0, expanding a(rj) in Taylor series 
around r] shows that a(rj) = 0. This means that the only singular bohmian trajectory is 
the trivial one of not having a universe at all! All non trivial quantum solutions have to be 
non singular. 



V. QUANTUM BEHAVIOUR OF A FLRW MODEL WITH DUST AND RADIA- 
TION 

As we have seen in section the superhamiltonian constraint for a FLRW model with 
non interacting dust and radiation is given by Eq. pijl : 

2 

W = -^--6Ka + ^+p v wO, (35) 
24a a v 

We see that both p v and p v appear linearly in TC, and their canonical coordinates r\ and ip 
are, respectively, conformal and cosmic time. As in the preceeding section, from the Dirac 
quantization procedure one obtains the quantum equation H^/ = 0, which reads 

[Wah ~ 6Ka ~ 'at, ~ %) * ( °' *> V) = °' (36) 
where we have used the usual coordinate representation p = —id/dq. 
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Either 77 or ip can be chosen as time parameters on which \1/ evolves. However, the classical 
solutions can be expressed explicitly only in conformal time rj [see Eq.(|29|)]. Furthermore, 
cosmic time ip depends on the constants characterizing each particular solution through 
tp — j drja(r]), and it is not the same parameter for all classical solutions (see Ref. 3] for 
deatils). Hence, we will take 77 (in fact —77, for the reasons mentioned in the previous section) 
as the time parameter of the quantum theory 1 . With this choice, and for a particular factor 
ordering, Eq. (j3T)j) can be written as: 



• 9 t, ^ 



(37) 



A. Eigenstates of total matter content 

In this subsection we only consider initial states |\l/(?7o)) which are eigenstates of the total 
dust matter operator p v . It follows that these states at time 77, \^(rj)) will also be eigenstates 
of p v with the same eigenvalue because [H,p^\ = 0. In other words, we consider that dust 
matter is not created nor destroyed. In such a way, we have p^l^ir])) = p v \ty(j])) and the 
wave function in the a, <p representation, (a,(p\^f(j])) = $?(a,(p,r)), is given by 



^From the Schrodinger's equation ([37)1. we have for ^>(a,rj) 



1 d 2 



+ Qua 2 — p^a \I/(a, rj), 



(38) 



(39) 



v 24 9a 2 j 
which is the Schrodinger equation for a particle of mass m = 12 in a one dimensional forced 
oscilator with frequency w = ^/K, and constant force p v , which we write as 



■ 9 T, \ 



1 d 2 



raw 



H —a -p v a ^{a,rj) 



(40) 



v 2m da 2 2 

The scale factor is defined only in the half line [0, 00), which means that the hamilto- 
nian (J35)) is not in general hermitian. Hence, if one requires unitary evolution, the Hilbert 
subspace is resctricted to functions on L 2 (0, 00; —00, 00) satisfying the condition: 



dip 



da 



a=0 



dip 



ggi (a,ip,T]) 
da 



*2 (a,<P,v) 



(41) 



a=0 



1 The choice of ip will probably yield a different theory, with a different Hilbert space. The kinetic term is 
more complicate and the measure is not the trivial one. We will not study this possibility here. 
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for any ^1(0, ip, rj), ^2(0, f, v) ^ L 2 (0, 00; — 00, 00). In the special case considered in this 
section, this condition is reduced to 



da 



(0 )77 ) =0^(0,77), 



(42) 



where a is a real parameter [25]. We will analyze the two extreme cases: a = and 
a = 00, which are the simpler and usually studied in the literature on quantum cosmology 
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16 



18, 



19 



21 



251. 



For the case a = we have that 



which is satisfied for even functions of a. For the case a = 00, we have 



*(0,7j)=0, 



(43) 



(44) 



which is satisfied for odd functions of a. Both of them address the boundary conditions of 
the wave packet near the singularity at a = 0. 

In order to develop the BdB interpretation, we substitute into the Schrodinger's equation 
fl4()j) . the wave function in the polar form \1/ = Ae tS , obtaining for the real part 



OS 1 (dS\ 2 mw 2 2 ^ n 



(45) 



where 



Q = 



1 d 2 A 



2mA da 2 

is the quantum potential. The Bohm guidance equation reads 



ma 



j dS_ 
da ' 



(46) 



(47) 



A solution ^(a, 77) of Eq. (j4()j) can be obtained from an initial wave function ^o(a) using 
the propagator of a forced harmonic oscillator. Let us do it for the two boundary conditions 
just presented. 
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1. The case of boundary condition a = 



We denote the propagator K a=0 (2, 1) = K a=0 (r]2, a 2 ; ai), where 1 stands for the initial 
time and initial scale factor rji, ai respectively, and 2 stands for their final values. 

The propagator when the Hilbert space is restricted to a > can be obtained from the 
usual one (i.e with coordinate — oo < a < oo) which is associated to a particle in a forced 
oscilator K(2, 1) = K(r] 2 , a 2 ; rji, a x ) by means of 

K a=0 {2, 1) = K( V2 , a 2 ; Vl , ai) + Kfo, a 2 ; Vl , -oi) (48) 

This symmetry condition is necessary to consistently eliminate the contribution of the nega- 
tive values of the scale factor j^J. The usual propagator associated to a particle in a forced 
oscilator is 



K(2,1) = J mW . exp(iS d ) (49) 
y 2m sm(wr]) 



where r\ = r\ 2 — r\\ . The classical action S c \ is given by 



Sd = tt^—, — A cos(wi])(al + a?) - 2a 2 a x + (a 2 + Oi) — ^[1 - cos(wr/)] - 
2 sin(w77j I " mw z 

[1 - cos(wry)] f y - + ^ y - sin(wr/)wr/l. (50) 
We assume that for 771 = 0, the initial wave function is given by 

*o(a) = (—) ' exp(-<xa 2 ), (51) 



7T 

where a > 0. The wave function in a future time r\ 2 is 

/•oo /-oo 

Ur(a 2) 77 2 )= / AT a=0 (2, ^^(aOda! = / tf(2,l)*o(ai)<fai, (52) 

JO ./-co 

where the even caracter of \l/(a, 0) has been taken into account to extend the integral. 
Integrating and renaming 77 = t] 2 , a = a 2 we have 



- a=0 / / VfT \ 1//4 / ( ""»' 



* a=u (a,77) = — J- ? r : T exp 



n J y i cos(wr])[2a tan(wri) — imw] l2tan(iu77) 



a 2 + 



mw ( p v \ 2 2ap v [1 - cos(wr?)] 

'-aH ^[1 - cos(w7/)J H — ? — + 



cos 2 (wr])[2a tan(wri) — imw] \ mw 2 J m w 2 cos(wr)) 

^pt ([cos(wrj) — 1] tan(w77)\ 



f[cos{wrj) — 1J tan{wr])\ 
m 2 V w 4 cos(wrj) ^ w 3 ) 



(53) 
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Flat spatial section (k = 0) 

We consider the case k — 0, which is obtained by taking the limit of the wave function 
given by Eq. (JHH)) for w — > 0. We compute its phase 5* from \1/ = Ae tS and calculate the 
derivative dS/da. In this way we have, for the Bohm guidance equation Eq. (j47)l . 

a' ^ a - 1 + „ (M) 

11 ~ A 2 2 i 2° ~~ //I 2 2 , T^PfV- l 04 J 



Comparing with the radiation case studied in [8J, we see that here it appears a term pro- 
portional to p v in the RHS of the Bohm equation. The general solution is: 



a{rj) = cJAg^ + w? + (55) 
v 2m 

where C is a positive integration constant. We can see that, contrary to the classical 
solution (}29|) . there is no singularity at r\ — 0. The quantum effects avoid it. Furthermore, 
for long times 7] 3> m/2a, Eq. ()55|) reproduces the classical behaviour (j2HJ) for the scale 
factor. 

For the case in which the evolution starts from a shifted gaussian wave function 

\&o( a ) = f~) ex P — CT ( a — a °) 2 ' (^) 
the Bohm guidance relation contains an additional term yielding the general solution 



0(77) = Co V4<rV + ™ 2 + fV + f , (57) 

which has exactly the same behaviour, apart from the shift on the minimal value of the scale 
factor by the ao/2 term. 

Positive curvature spatial section (k — 1) 

Setting w = = 1 in the wave function given by Eq. (|B^j) and computing its phase S, 
we obtain for the bohmian trajectories 

a(rf) = Cq\J Aa 2 sm 2 (i]) + m? cos 2 (r]) + 7r~\\ ~ cos ( 7 ?)]; (58) 

where Co is a positive integration constant. This is a non singular cyclic universe, see 
figure 2, which presents classical behaviour for rj such that | tan(7/) |3> m/2 [see Eq. (}2T?j) ]. 
Quantum effects avoid the classical big bang and big crunch. 
Negative curvature spatial section (k = —1) 
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FIG. 2: Scale factor evolution of a quantum universe with positive curvature spatial section filled 
with matter and radiation. 

Setting now w = = i in the wave function yields the bohmian trajectories: 



a(r)) = Cq\J Aa 2 sinh 2 (?7) + m 2 cosh 2 (77) + -^-[cosh(?7) — 1]. 



(59) 



Again, Cq is a positive integration constant. This is a non singular ever expanding universe 
which presents classical behaviour for 77 such that | tanh(ry) |3> m/2 [see Eq. ()29j) ]. Quantum 
effects avoid the classical big bang. 

As in the k = shift in the center of the initial gaussian will not modify these 

solutions qualitatively. 

For the boundary conditions a = 00, or \l/(0,t) = 0, the propagator K a= °°(2, 1) can be 
obtained from the usual (i.e, with coordinate —00 < a < 00) propagator associated to a 
particle in a forced oscilator K(2, 1) by means of 



In order to satisfiy the condition \I/ (0, 77) = 0, we take as the initial wave function a wave 
packet given by 



where a > 0. Following a similar procedure as in the case a = 0, we calculate the wave 
function by propagating the initial wave function as 



K a -°°(2, 1) = K(t] 2 , a 2 - 771, ai) - K(r) 2 , a 2 ; 771, -a x ). 



(60) 




(61) 
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poo poo 

y(a 2 , V2 )= K a=oo (2,l)^ (a 1 )da 1 = K(2,l)S> ( ai )da u (62) 

JO J-oo 

where the odd caracter of ^ has been used in order to extend the integral. Integrating this 
expression and renaming a = a 2 and r] = r]2 with rji = 0, we have 

^=^(a,r ] ) = (^)^=°(a, V ) (63) 



where 



_ vmw 



sm(wr]) 
and 



—a H — ^t(1 — cos(wr]) 
mw 2 



(64) 



L> = 7 — r - a 65 

2tan(wr/) v ; 

The phase of \1/ Q=00 (a, 77) can be expressed as the sum : 

phase[^ a= °°(a,r/)] = P hase (^") + phase[^ a=0 (a, 77)], (66) 

and it is easy to see that the phase of (-C/2D) is independent of a. Then, 
[<9phase(\I/ Q=00 (a, rj))]/da = [<9phase[\I/ a=0 (a, r])]]/da, and the Bohm guidance relations are 
the same as in the previous cases. Therefore, the solutions are the same. 

The quantum cosmological models obtained in this subsection have the nice properties of 
being non singular and presenting classical behaviour for large a. However, they suffer from 
a fundamental problem: the wave function (J53|) from which they are obtained does not have 
an unitary evolution. The reason is that propagators constructed from Eq.'s ()48j) and ()6()|) 
do not in general preserve the hermiticity condition (J42j) imposed on the wave functions: 
it depends on the classical potential. In Ref. Q], there are obtained the potentials which 
allow propagators in the half line (a > 0) to preserve unitary evolution. The potentials of 
the previous section are some of them but the potentials of the present one are not. Hence, 
even though the initial wave function Eq. (J51|) satisfies the hermiticity condition, the wave 
function (|53|) does not. Let us then explore the more general case of initial superpositions 
of the total dust mass operator eigenstates. 
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B. Analysis of wave packets given by superpositions of total dust mass eigenstates 



In this subsection we consider the case of a general solution of Eq.(|37|) which is not 
necessarily one of the eigenstates of p v , the total dust mass operator. 

Following the BdB interpretation of quantum mechanics, we substitute in Eq. (J37|) the 
wave function in polar form: \1/ = A (a, <p, 77) exp {iS (a, tp, 7])}. The dynamical equation 
splits in two real coupled equation for the two real functions S and A (recall that w = ^[k 
and m = 12). 

dS 1 fdS\ 2 dS mw 2 2 ^ n 

^ + 2^(*j -% + — ° +<? = °- (67) 

^d! + A ( a a 2 ) + — (a 2 ——) - (68) 
dr] d^p ^ ' da \ m da J 

where 

Q s -tkjiw- (69) 

Equation (|67p is the modified Hamilton- Jacobi equation where Q (a, <p, 7]) is the quantum 
potential which is responsible for all the peculiar non classical behaviours. When the quan- 
tum potential is zero, the equation is exactly the classical Hamilton- Jacobi equation. The 
momenta are given by the Bohm's guidance equations 

A - (70) 

Note also that 

» - ™^ (72) 

is the total 'energy' of the system, which is interpreted, from its classical meaning, as the 
total amount of radiation in the universe model. 

In the causal interpretation, equation (|HH)l is a continuity equation where A 2 is a proba- 
bility density. The generalised velocities can easily be identified as 

a' es lgg%M, (73) 
m da 

if' = a. (74) 

Consider now the classical limit (Q = 0). Then the solution of the principal Hamilton 
function (S) is just S = W (a) — Ei] + p^ip, where E and p v are constants. Since p v is 
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da 


dS 


(a,ip,r]) 




dip 


dS 


{a,<p,ri) 



proportional to the total amount of dust matter in the universe, and E to the total amount 
of radiation, there is no creation or annihilation of dust matter nor radiation. However, in 
the presence of a quantum potential, this solution is no longer valid, opening the possibility 
of non conservation of matter and radiation due to quantum effects. 

1. Formal Solutions 

We now turn to the problem of solving the Schrodinger's equation (|37|). 
Flat spatial section 

For the case of flat spatial section (k = 0), equation (pTTjl simplify to 



.d^f (a,(p,rj) 1 d 2 ty(a,ip,r]) . (a,(p,r)) 



drj 2m da 2 

To solve this equation we make the ansatz 



dip 



* (a, <P,v) = X (a) exp [~^P v) ex P *A 
where x ( a ) must satisfy the differential equation 2 



(75) 



(76) 



d 2 X (a) 
da 2 



+ vax (a) + fix (a) = 0. 



This is essentially an Airy equation with solution given by 



(77) 







la + - <AZi 
v I 3 



BZ 



(78) 



The Zi function is the first kind Bessel function of degree ~, and the A and B can be any 
functions of v and (3. 

The general solution is a superposition given by 



x <A(p,v) Zi 



+ B(P,v) Z_i 



2 As in the following we will make superpositions of eigenfunctions of the total dust matter operator, we 
will use from now on the letter v in order to not confuse it with the beable p v — dS/dip. We did not 
make this distinction before because they coincide for eigenfunctions of the total dust matter operator. 
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Positive curvature spatial section: Landau levels 
In the positive curvature case (k—1), Eq.(pT7|) reads 



- + — aw (a, (p, rj) + ta 



<9w (a, ^,77) 



(79) 



9?7 2m 9a 2 2 " <9<£> 

There is a canonical transformation which simplifies the problem. Let us define new variables 
given by 



£ = \pma 



P<p 



m 



a = —y/mip + 



Pa 



Pa 

Pi = -?= ; p° 



p<p 



m 



m 



Using these new variables, the hamiltonian decouples in two parts, one describing a harmonic 
oscillator and the other a free particle: 



H 



5(* + ? 



1 



(80) 



harmonic oscillator free particle 
Decomposing the wave function as 



*(t,o-, V )=x(0exp{-i(er ] + V2ka)}, (81) 

we immediately recognize that x (0 = ex P { — } hn (£) , where h n are the Hermite polino- 
mial of degree n. Just as for the harmonic oscillator, the index e is constrained to take the 
values 

e n = k + [n + 



(82) 

where k can take any real positive value while n is apositive integer. Eq. ()82j) determines 



a set of Landau levels for the cosmological model 
superposition given by 



22|. The most general solution is a 



oo . 

*(£,<t,77) = E /dA:xn(e)[^n(fe)exp{i < 7V^fe} + 

n=CT 

G n (k) exp |— icr \/2A;} x exp {— i e n 77} . 



(83) 



The quantities D n (k) and G n (k) are arbritary coefficients that can depend on the param- 
eter k. Recall that we have performed a canonical transformation that mix coordinates 
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and momenta, and these are not the proper variables to apply the causal interpretation. 
Instead, it is imperatif to apply the inverse transformation to the coordinate basis before 
using the guidance relations. This is a necessary requirement to maintain the consistency 
of the causal interpretation of quantum mechanics 27 1- 2sj]. 

For the negative curvature spatial section (re = — 1), the general solutions are hypergeo- 
metric functions whose asymptotic behaviours are rather complicated to study in order to 
obtain reasonable boundary conditions. Hence, we will not treat this case here. We proceed 
to the analysis of an interesting particular solution. 



2. Transition from exotic dust to dust in the flat case 



The quantum states of the matter and radiation FLRW universe studied in section IV A II 
are eigenstates of the total dust matter operator p v . The total wave function is given by 
\I/(a, (p, rj) = \l/(a, rj) exp(iipp 9 ) where \I/(a, rj) is given by Eq. (Jo"3|) . Taking the limit w — > in 
that equation, we obtain the wave function ^/(a,r]) for the case of flat spatial section, k = 
which, after renaming the eigenvalues of total mass by v = p v , is given by 



(a, rj) 



Sam 



7T fx 



exp 



m 
+i — 
2r] 



where 



m 2 a 



/i 



v rj 
2m 



2m 



a 2 2 i 2 

jj, = 4a rj + in , 

a = arctan . 

\ m J 



v rf" 6 

i i — h 

6m 2 



vrf\ 2 m 2 ( vt] 
a + — — I I a 



2m 



(84) 



Now we consider a more general situation than in section IV A II We suppose an initial 
state at rj = which is given by a gaussian superposition of eigenstates of total matter 



/OO 2 
dv exp- 7( ^ o) *„(a,0) exp{-z>i>}. 
-oo 

Then, the state at time 77 is given by 

dv exp~' 1[v ~ V0 > $! v (a,r}) exp{-i(pv}. 

-00 



(85) 



(86) 



21 



In this way, we have a square-integrable wave function. We find 



■5a nm 



i 

2\ 4 



exp 



'ft(F) am 2 \ 2 3? (J) 2 3fe(L) 



4i/ 



a + 



where we defined 



3? (M) 3? (P) 



4/y 



4z/ 



'3(F) 
4z/ 



+ 



4i/ 
m 

2/J,T] 



a if 



Av 



Av 



a+ 



{[i - m 2 ) a 2 H j-^a v?+ 



4z/ 



4// 



4z/ 



4i/ 



7 + 



err/ 



+ 



r = arctan 



4/i / (24m//)' 
r/ 3 (/i + 3m 2 ) 
24m(7/i + err/ 4 ) 



(/i + 3m' 



F 
G 

J: 



mail ii / 9 

/i + m 2 

/i 2/i v 



ail . r/ 

7+ — ' 



4/i 24m/i 



f/i + 3m' 



-2z 



mart n / 

L + + 

/x 2/i v 



rrr 



7 + 



r/° 



err/ 



4/i 24m/i 



( // + 3m' 



7 + 



err/ 



4/i 24m/i 
L = —2i 7?j G 

M = 4 r 7?Jo J 
p = -4 7 2 tj 2 J 



ux + 3m' 



and 9? and 3 stands for the real and imaginary part, respectively. 

If one calculates the squared norm of the wave function, one obtains 

'da I" dip\m\ 2 = J— l + 4=erff^!' 

/0 J-oo V 7 \/7T V 2m 



/"CO /"CO 

/ da df\\^\ 

JO J-oo 
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(87) 



where erf (a;) is the error function. The only dependence on time can be eliminated by 
choosing the gaussian to be centered at vq = 0. With this choice we garantee unitary 
evolution of the total wave function. ^From equations (|71| ) - (j74j) . the trajectories can be 
computed by solving the given system of equations 



2 

m 



3(F) 
Av 



+ 



m 



2/ir/ 



(u — m 2 ) 



a | *M 
Amu ^ 



(88) 
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ip = a 



(89) 




3(G) 



a 



(90) 



Note that p^ is no longer constant. We integrated numerically these equations with the 



close to the origin as shown on figure 3. Hence, we expect that quantum effects be relevant 
only in this region. Far form the origin, the scale factor must behave classically. 



FIG. 3: Quantum potential as a function of time. For regions far from the origin, the quantum 
potential goes to zero, implying a classical behaviour. 

The behaviour of p v is plotted on figure 4. ^From this plot we can see that far from the 
origin p v is constant. This is in accordance with classical behaviour as long as the quantum 
potential is zero in this region. The surprising feature is that in the far past the universe 
was filled with a classical exotic dust (p v < 0). 

^From Eq. ()72j1 . one can also compute the amount of radiation. Figure 5 shows the result. 
Again, far from the origin, radiation is conserved while in the origin, due to quantum effects, 
it is not conserved. 

For the evolution of the scale factor, numerical integration of equation ()88|) yields the 
plot of figure 6. In the far positive region, the scale factor behaves classically, as expected, 
and matter is conserved. On the other hand, in the far negative region the scale factor 
also behaves classically but with a universe filled with exotic dust, and here again matter 
is conserved (compare this region with figure 1). Both regions have a consistent classical 
behaviour. Hence, the universe begins classically from a big bang filled with exotic dust 



renormalisation condition a (0) = 1. The quantum potential Q 



1 8 2 A 
2m A da 2 



is non zero only 



.Quantum Potencial 
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FIG. 4: Evolution of the total amount of matter in the universe (p^). In the far past, the uni- 
verse was filled with exotic dust (p^ < 0). Close to the origin, quantum effects transform it into 
conventional dust (p^ > 0). 



FIG. 5: Evolution of the total amount of radiation in the universe (p v )- Far from the origin the 
universe behaves classicaly and the total amount of radiation is conserved. It varies near the origin 
due to quantum effects. 

and conventional radiation. It evolves until it reaches a configuration when quantum effects 
avoid the classical big crunch while transforming exotic dust into normal dust. From this 
point on, the universe expands classically filled with conventional dust and radiation. 




Total Amount of Radiation 

5- 



-10 



-8 



10 



24 




FIG. 6: Evolution of the scale factor. Close the origin quantum effects avoids the collapse. For the 
others regions, the evolution is essentially classical. 

VI. CONCLUSIONS 

In the present work we studied some features of the minisuperspace quantization of FLRW 
universes with one and two fluids. For the one fluid case (radiation), we have generalized 
results in the literature by showing that all bohmian trajectories coming from reasonable 
general solutions of the wave equation obtained through the assumptions of unitarity and 
analyticity at the origin, do not present any singularity. Hence, this quantum minisuperspace 
theory is free of singularities. 

For the two fluids case (non interacting radiation and dust), we first obtained bohmian 
quantum universes free of singularities reaching the classical limit for large scale factors. 
However, these trajectories arise from eigenfunctions of the total dust mass operator whose 
time evolution is not unitary. When considering the general case, we managed to obtain 
a wave solution presenting unitary evolution with some surprising effects. Now dust and 
radiation can be created but the new feature is the possibility of creation of exotic fluids. 
We have shown that dust matter can be created as a quantum effect in such a way that 
the universe can undergo a transition from an exotic dust matter era to a conventional dust 
matter one. In this transition, one can see from figure 5 that radiation also becomes exotic 
due to quantum effects, helping the formation of the bounce. 

The fluid approach is not fundamental, but we expect that it can be quite accurate in 
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describing quantum aspects of the Universe, in the same way the Landau description of 
superfluids in terms of fluid quantization was capable of showing many quantum features of 
this system j^. After all, creation and annihilation of particles as well as quantum states 
with negative energy are usual in quantum field theory. The formalism developed in the 
present paper seems to be a simple and calculable way to grasp these features of quantum 
field theory. Their physical applications may be important: exotic fluids are relevant not 
only in causing cosmological bounces and avoiding cosmological singularities jsij], but also 
for the formation of wormholes 0, and for superluminal travels j^]. These are some 
developments of the present paper we want to explore in future works. 
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